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I. Introduction 

The asymptotic behavior of adaptive systems has been a well researched topic during the last 
couple of decades, and it is well known that asymptotic tracking can be achieved using the Lyapunov 
redesign method. However, the transient behavior of the input and output signals can be very 
oscillatory with significant excursions [1], There has been a great deal of effort to modify the 
control architecture and the adaptive laws from the perspective of improving the transient behavior 
of the tracking error. The majority of these efforts led to nonadaptive high gain feedback [2, 3], 
switching control law [4] or to a parameter dependent persistent excitation condition [5]. 

On the other hand, the tracking error magnitude can be decreased by increasing the adaptation 
rate, but this results in unwanted high frequency oscillations and large overshoot in the control signal. 
This shortcoming is common for the majority of existing adaptive control methods. Recently some 
results have been appearing in the control community, which explicitly address the input signal 
transient behavior (see for example [6]). 

In this note, instead of modification of the control architecture or the adaptive laws, the reference 
model is modified by feeding back the tracking error signal. This approach, called here "modified- 
reference- model MRAC" or M-MRAC in short, is motivated by the fact that the initial large errors in 
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the control gains generates large transient excursions both in the system’s input and output signals. 
The idea is to prevent the system’s attempt to aggressively maneuver toward the reference model 
by modifying the reference model by a term proportional to the tracking error. As the tracking 
error approaches zero, the reference model approaches its original form. Therefore, the system 
asymptotically tracks not only the modified reference model, but also the original one. In addition, 
the error feedback term determines the damping in the control signal, which in concert with the 
adaptation rate makes it possible to regulate the transient of the control signal. A design guideline 
is provided for the selection of the feedback gain relative to the adaptation rate. The proposed 
adaptive control method has uniform performance in reference commands and initial conditions 
without the need for re-tuning. 

It is worth to note that the reference model modification is not new. It was used to accommodate 
for the input saturation in the adaptive system (see for example [7, 8] and references therein). 

II. Reference model modification 

Consider a multi-input multi-output uncertain linear system 

x(t) = Ax(t) + Bu(t), x(0) = xo , (1) 

where x £ 1” and u £ R 9 are the state and input of the system respectively, and A £ R nxn and 
B £ M. nxg are unknown constant matrices satisfying the following matching conditions. 

Assumption II. 1 Given a Hurwitz matrix A m £ M” x " and a matrix B m £ R nX9 of full column 
rank, there exists a matrix K\ £ R raX9 and a sign definite matrix A £ R <JXlJ such that the following 
equations hold 

B = B m A, A = A rn — BKj . 

Remark II. 1 The sign definiteness of A corresponds to the conventional condition on the high 
frequency gain matrix of MIMO systems (see for example [9]). Without loss of generality it is 
assumed that A is positive definite. The matching equations (2) complete the conditions for the 
existence of an adaptive controller. 
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As in the conventional model reference adaptive control (MRAC) framework, the objective is 
to design a control signal u(t) such that the the system ( 1 ) tracks a given reference model 

x r (t) = A m x r (t) + B m r(t), x r (0) = xq , 

where A rn . B m are chosen according to performance specifications and satisfy Assumption II. 1 and 
r(t) is a bounded and piecewise uniformly continuous external command. It can be noticed that 
discontinuous commands such as steps or square waves satisfy these conditions. 

The system (1) can be written in the form 

x(t.) = A m x(t) + B m r(t) + B m A[u(t) - 0 T /(t)] , (2) 

where 0 T = [A'-, 1 A~ T ], f(t) = [£c T (t) r T (t)] T . Ideally, if 0 were known, one could set the control 
signal equal to 

u°(t) = 0 T /(f) , (3) 

thus reducing the system’s dynamics to the given reference model (2). Although the reference model 
(2) can always be specified from the performance perspectives, the control signal (3), which here is 
called a reference control signal, can not be implemented. Therefore, the adaptive version of it is 
implemented in the MRAC framework, that is 

= (4) 

where 0 (f) is the estimate of the ideal control gain 0 , and is updated online according to the 

adaptive law 

6 (f) = -7 f(t)e T (t)PB m , 0(0) = 0 O , (5) 

where 7 > 0 is the adaptation rate, e(t) = x(t) — x r (t) is the tracking error, and P = P T > 0 is 
the solution of the Lyapunov equation 

Aj n P + PA rn = -Q (6) 

for some Q = Q T > 0. Introducing the parameter estimation error as 0(f) = 0(f) — 0, the tracking 
error dynamics can be written in the form 

e(t) = A m e(t) + R m A0 T (f)/(t) . (7) 
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It is well known that this control architecture guarantees asymptotic tracking x(t) —> x°(f) as 
t, — > oo, while ensuring the boundedness of all closed- loop signals. However, the transient behavior 
of x(t) and u(t) cannot be guaranteed, because high frequency oscillations are generated in the 
control signal, when the adaptation rate is increased in order to obtain better tracking in transient. 
The proposed M-MRAC architecture employs a modified reference model in the form 

x m (t) = A m x m (t ) + B m Nr(t) + A e m (t), x m (0) = x 0 , (8) 

where A > 0 is a design parameter and e m (t) = x(t) — x m (t) is the modified tracking error, and the 
subscript in is introduced to distinguish between the original reference model (2) and the modified 
reference model (8), which is used for the actual control design. The original reference model and 
the reference control signal will be used only for the analysis. The modified tracking error satisfies 
the dynamic equation 

e m {t) = ( A m - AI„)e m (£) + B m A[u(t) - 0 T /(f)] , (9) 

where I„ denotes ?r-dimensional identity matrix. The adaptive control u(t) is defined according to 
equation (4) as in MRAC design, but the adaptive law now is based on the modified tracking error 

^m(£) 

kt) = -7 f(t)el(t)PB m , 0(0) = 0o . (10) 

With the application of the adaptive control u(t) the modified tracking error dynamics reduce to 

e m (t) = (A m - AI n)e m {t) + B m A@ T (t)f{t) , (ff) 

where the parameter estimation error 0(f) is defined similar to MRAC. 

In the following analysis, the control error signal is needed, which is defined as u(t) = u(t) — 
u°(t). From the above constructions it follows that 

u{t.) = 0 T (t)/(t) . (12) 

Since the reference control signal is the best achievable signal, it is in our interest to minimize the 
control error u(t), as well as the state errors e(t) and e m (£), both in transient and steady state 
by selecting proper values for the adaptation rate 7 and feedback parameter A. This is the main 
objective of the analysis in the following sections. 
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Remark II. 2 It is worth to point out that the modified reference model is used for the control design 
purposes. The asymptotic and thransient analysis of the adaptive system is conducted with respect 
to the original reference model, which implies that the desired dynamic characteristics built in the 
specified reference model are preserved with the proposed control method. 

III. Analysis of M-MRAC Design 

A. Asymptotic properties 

The following theorem summarizes the asymptotic properties of the proposed M-MRAC archi- 
tecture. 

Theorem III.l Let the system (1) satisfy Assumption II. 1, and the reference model be given by 
(2). Let the system be controlled by the M-MRAC scheme given by (f), (8) and (10). Then 

1 ) all closed-loop signals are bounded 

2) e m (t) — > 0, eft) — > 0 and uft) —>0 as t — > oo. 

The theorem can be proved using the standard Lyapunov tools and the extended Barabalat’s 
lemma from [10]. 

Theorem III.l shows that M-MRAC architecture guarantees asymptotic tracking of not only 
the modified reference model (8), but also the original reference model of the conventional MRAC 
architecture (2). That is the asymptotic performances of both conventional MRAC and proposed 
M-MRAC designs are equivalent. 


B. State transient properties 

Using the Lyapunov analysis it can be shown 


I -Coo < 


\Arni „(P) Vl ’ 


where a = trace ^A0 T (O)0(O)^ and 


I e m (Oil £2 — 


\J Amin (Q ) + 2AA min (P) \/7 
The bounds on the tracking error signal eft) are obtained using the equation 


(13) 


(14) 


— [e(t) - e m ft)] = A m [e(t) - e m ft)] + A e m ft) 


(15) 


5 




C. Input transient properties 

It can be noticed that u(t) does not explicitly depend on design parameters A and 7. Instead, 
u(t) depend on 7 through the adaptive laws, and u(t) depend on A through the tracking error 
dynamics. Therefore we conduct dynamic analysis of this signal. It satisfies the differential equation 

u{t) + A u(t) + 7 p(t)Gu(t) = 77/ {t)B^Pe m {t) + (s + A )h(t) , ( 17 ) 

where p(t) = G = B^PB m A, rj(t) = -p(t)Bj l PA m e m (t) - p{t)B^Pe m {i:) and h{t) = 

Q T (t)f(t). Since all closed-loop signals are bounded, it follows that ||p(t)||,c 00 < ai, ||p(t)||£ 00 < «2 
and ||/i’(t)||£ ao < a 3 for some positive ai, «2 and a.3. Therefore, the equation ( 17 ) can be considered 
as a second order linear equation with time varying coefficients in u(t). It can be inferred from the 
equation ( 17 ) that the adaptation rate 7 determines the frequency of u(t) and hence the frequency 
of the control signal u(t), since the reference control u°(t) is in the low frequency range. Therefore, 
increasing 7 increases the oscillations in u(t) as it is the case for the conventional MRAC design. On 
the other hand, A determines the damping ratio. Therefore increasing A suppresses the oscillations 
in u(t) and hence in the control signal u(t). 
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A proper value for the parameter A is selected from the perspective of minimizing the 
norm bound on u(t). Following the steps from [11], it can be shown that selection of A = 
y / 2ai7A max (-Bm-PI? m ) results in the bound 

||«(*)|| < ke~ ut +/3 4 7“® , (18) 

where the positive constants /?3 and p 4 do not depend on 7 and v is proportional to ^7. 

Using the inequality (18), a tighter bound can be obtained for tracking error form the error 
dynamics (7). This bound has a form of 

||e m (t)|| < I3 5 e~ , '* t +/3e7 _ ® , (19) 

where v* = rninfzz, u m ] and v m is the rate of decay of e Amt . 

It follows from the inequalities (18) and (??) that the input and output transient errors can be 
decreased as desired by increasing the adaptation rate. 

IV. Application to a Generic Transport Model 

In this section, the advantages of the M-MRAC architecture are demonstrated in simulations for 
a dynamic model that represents the lateral-directional motion of a generic transport aircraft (GTM) 
[12]. The nominal model is the linearized lateral-directional dynamics of GTM at the altitude of 
9144 to and speed of 0.8 M and is given by the equation 

x(t) = A n x(t ) + B n u(t) , (20) 

where x = [/? r p </>] T is the lateral-directional state vector, in which /3 is the sideslip angle, r 
is the yaw rate, p is the roll rate, </> is the bank angle, and u = [d 0 d r ] T is the control signal that 
includes the aileron deflection 6 a and the rudder deflection S r . The controlled output is selected to 
be the sideslip angle and the bank angle of GTM. The open loop system has a slow Roll mode and 
negligible damping in the Dutch- Roll mode (see Table 1 for details). 

The reference model is selected from the perspective of improving the performance character- 
istics of the nominal dynamics and is given by the matrices A rn = A n — B n Ko and B rn = B n N, 
where the feedback and feedforward matrices are selected to speeds up the Roll mode and increases 
the damping ratio of the Dutch-Roll mode. 
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Table 1 Dynamics characteristics of the nominal model. 


Open-loop Eigenvalues 

-0.0177, -1.6339, -0.3212 ± 1.7404i 

Closed-loop Eigenvalues 

-1.1623, -3.5583, -1.2106 ± 1.2008« 


The uncertain model of GTM corresponds to about 28% loss of left wing tip and 55% loss of 
rudder surface. The control objective is to track a reference command r(t) = \/3 CO m(t) (f>com(t)] T 
form the zero initial conditions. The reference command is chosen to be a series of coordinated turn 
maneuvers. That is, /3 CO m{t) = 0 and the bank angle command is chosen to be a square wave of the 
amplitude of 15 degrees and of the frequency jjj rad/sec. In order to make (j> CO m(t ) differentiable, it 
is filtered through a first order stable filter s + 1Q . For each 7, parameter A is selected with oq = 

First, a simulation is run with 7 = 100 and Q = I4. Figure 1 displays the output tracking 
of the adaptive system with M-MRAC and MRAC versus the output of the original reference 
model. Clearly M-MRAC has a superior performance even for the chosen low adaptation rate. The 
corresponding adaptive control signals and the reference control signals are displayed in Figure 2. 
It can be observed that the M-MRAC control signals tracks the reference control signal quite well 
without any oscillations unlike the MRAC control signal. 



Fig. 1 Output performance of the proposed M-MRAC with 7 = 100. 

Second, it is demonstrated that both the output and input tracking errors can be decreased by 
increasing 7, without generating the oscillations in the control signal. For this purpose, simulations 
are run with 7 = 40000 and 7 = 200000. The results are presented in Table 2. As it can be seen 
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Fig. 2 Input performance of the proposed M-MRAC with 7 = 100. 

from the table, the control errors are systematically decreased when 7 is increased in the case of 
M-MRAC design. In the MRAC design case these errors are large for large values of 7. The tracking 
errors are decreased when 7 is increased for both designs. In general, the tracking error may be 
smaller in the case of MRAC, since the term f ,c o J on the right hand side of the first inequality 
in (16) is absent in the case of MRAC design (cq=0). 


Table 2 Output and input tracking error norms in degrees 
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Design 


(l) 1 1 Cqq 

Sa(t) 

Sr(t) 

400 

MRAC 

0.2826 

0.7513 

15.0577 

23.0739 

M-MRAC 

0.0255 

1.3035 

11.7093 

0.6349 

40000 

MRAC 

0.0056 

0.0172 

27.2261 

47.1725 

M-MRAC 

0.0036 

0.1958 

4.6656 

0.2307 

200000 

MRAC 

0.0025 

0.0051 

27.3205 

47.1719 

M-MRAC 

0.0016 

0.1162 

3.7371 

0.1922 


Next, to demonstrate the scaling properties of M-MRAC, a simulation is run with 2-times scaled 
< t>com{t )• In this case, the output performance does not change, as can be seen from Figure 3, and 
the input is scaled 2-times as in the case of a linear time invariant system. 

The final simulation is run with nonzero initial conditions to show that M-MRAC does not 
degrade with the change of initial conditions, unlike conventional MRAC. The output and input 
performance of M-MRAC with 7 = 200000 and xq = [0.2 0.5 — 3.0 10] T is displayed in Figures 
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Fig. 3 Output performance of the proposed M-MRAC for scaled command. 
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Fig. 4 Input performance of the proposed M-MRAC for scaled command. 

5 and 6. 

It can be observed that a good tracking of the output and input signals is achieved both in 
transient and in steady state. 


V. Conclusions 

A simple yet powerful modification of the conventional MRAC design is presented, which has 
a guaranteed transient and asymptotic performances in both the input and output signals of the 
system. The proposed M-MRAC design enables one to achieve close tracking of reference output 
and input signals by increasing the adaptation rate without generating high frequency oscillations 
in the control signal. The performance of M-MRAC is uniform in initial conditions and in reference 
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Sideslip angle in degrees 



Bank angle in degrees 



Fig. 5 Output performance of the proposed M-MRAC for nonzero initial conditions. 



Fig. 6 Input performance of the proposed M-MRAC for nonzero initial conditions. 

command. 

Since only the reference model is modified leaving the adaptive laws intact, this method can 
be applied along with any known adaptive laws, such as dead-zone modification, e-modification, 
cr-modification, etc. 
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